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Scripps  Institution  of  Oceanography 
University  of  California,  San  Diego 
San  Diego,  CA.  92152 


ABSTRACT 

When  explosive  charges  are  used  as  sound  sources  in  marine  refraction  experiments,  a  streamer 
hydrophone  is  often  trailed  near  the  surface  by  the  vessel  which  launches  the  charges.  The  output  of  the 
streamer  hydrophone  is  used  to  monitor  the  detonations,  and  after  correcting  for  the  time  spent  in  traveling 
from  the  shot  to  the  streamer  hydrophone,  provides  a  time  origin  for  each  shot.  If  three  arrivals  are  present 
(the  direct  water  path  and  two  reflections),  then  information  about  the  shot  depth  and  sea  floor  depth  may 
also  be  obtained. 

Since  the  ship  that  launches  the  charges  generally  steams  away  from  the  drop  site,  incorrect  results 
are  obtained  if  vertical  raypaths  are  assumed.  In  this  report,  a  method  which  does  not  assume  vertical  ray- 
paths  is  derived  to  find  shot  and  sea  floor  depths.  Inputs  are  the  distance  traveled  from  the  drop  site  when  the 
detonation  first  appears  on  the  streamer  hydrophone  output,  and  the  time  differences  beiween  uie  rnree 
arrivals.  A  Fortran-77  program  which  implements  the  method  is  listed  and  its  use  is  illr.  rr.:::'.  on  data 
obtained  during  a  refraction  experiment  conducted  on  the  Monterey  Deep-Sea  Fan. 
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Iterative  method  for  the  estimation  of  shot  and  sea  floor  depths  using 
hydrophone  streamer  data 


1.  Introduction 

Marine  refraction  experiments  are  conducted  using  an  impulsive  airgun  or  explosive  charges  as 
a  sound  source.  A  typical  experiment  (Figure  1)  might  have  one  vessel  utilized  as  a  receiving  station 
with  a  hydrophone,  or  a  hydrophone  array,  deployed  in  the  water  column.  A  second  vessel  (called  the 
shooting  ship)  steams  away  from  the  receiving  ship  and  drops  explosive  charges  into  the  water.  Acoustic 
energy  travels  along  numerous  paths  in  going  from  the  shot  to  the  array.  Since  there  is  a  general  increase 
in  sound  speed  with  depth  in  the  sediment,  much  of  the  energy  that  enters  the  sediment  is  refracted  back 
into  the  water  column  and  received  at  the  array.  Analysis  of  this  received  data  yields  information  about 
the  velocity  and  attenuation  structure  in  the  sediment. 

When  explosive  charges  are  used,  a  hydrophone  (called  a  streamer  hydrophone)  is  often  trailed 
near  the  surface  by  the  shooting  ship  while  charges  are  being  launched.  Streamer  hydrophones  are  used 
during  marine  refraction  experiments  to  monitor  the  detonations.  Output  from  the  streamer  hydrophone 
is  telemetered  to  the  receiving  ship  via  a  radio  voice  channel  and  recorded  along  with  data  from  the  array. 
Streamer  hydrophone  output  indicates  when  the  shot  detonates,  and  after  correcting  for  the  time  spent  in 
traveling  from  the  shot  to  the  streamer  hydrophone,  provides  a  time  origin  for  each  shot.  Additionally,  if 
three  arrivals  can  be  seen  on  the  streamer  hydrophone  output,  estimates  of  the  shot  and  sea  floor  depths 
may  be  obtained. 

An  example  of  this  type  of  data  is  given  in  Figure  2.  The  top  nine  time  series  are  the  array  data 
and  the  bottom  time  series  is  the  output  of  the  streamer  hydrophone.  In  the  streamer  data,  the  pulse  at 
0.32  seconds  is  the  arrival  which  travels  from  the  shot  directly  to  the  streamer  hydrophone.  The  second 
arrival  (at  1.51  seconds)  is  a  bottom  reflection,  and  the  third  arrival  (at  3.69  seconds)  is  a  surface 
reflection  -  bottom  reflection.  Energy  appears  at  the  array  approximately  14.7  seconds  after  it  arrives  at 
the  streamer  hydrophone. 

Ray  paths  corresponding  to  the  arrivals  on  the  streamer  channel  are  illustrated  in  Figure  3(a). 
The  first  arrival  is  denoted  (D),  the  second  arrival  (B),  and  the  third  (SB).  The  shooting  ship  moves  a 
horizontal  distance  x  from  the  drop  site.  The  water  column  is  modeled  as  a  two  layer  medium.  Layer 
boundaries  are  specified  by  the  surface,  shot  depth  ( d ])  ,  and  sea  floor  depth  (d\  +  d^-  Each  layer  is 
assumed  to  have  a  constant  sound  speed  (v[  and  vj).  This  is  not  a  very  good  model  of  a  velocity  profile 
for  the  water  column,  but  if  mean  values  for  each  layer  are  found  from  the  velocity  profile,  the  method 
gives  satisfactory  results. 
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Figure  1.  Refraction  experiment. 

A  hydrophone  array  is  deployed  at  mid-depth  in  the  water  column.  The  shooting  ship  steams 
away  from  the  array  and  drops  explosive  charges  into  the  water.  There  is  a  general  increase  in 
sound  speed  with  depth  in  the  sediment,  so  much  of  the  acoustic  energy  that  enters  the  sediment 
is  refracted  back  into  die  water  column.  Analysis  of  this  data  yields  information  about  velocity 
and  attenuation  in  the  sediment. 
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Figure  2.  Hydrophone  array  and  streamer  data. 

The  top  9  time  series  are  array  data,  bottom  time  series  is  the  output  of  the  streamer  hvdrophone 
Three  arrivals  are  present  on  the  streamer  channel  (at  .32,  1.51,  and  3.69  seconds).  Energy  appears 
at  the  array  approximately  14.7  seconds  after  the  First  arrival  at  the  streamer  hydrophone 
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2.  Derivation 


If  the  angles  in  Figure  3(a)  are  small  (i.e.  the  ship  does  not  travel  very  far  from  the  drop  site), 
then  the  raypaths  are  almost  vertical,  and  the  time  difference  between  the  second  and  third  arrivals  (At  2,3) 
is  due  to  the  extra  distance  that  the  (SB)  ray  travels  in  layer  1.  In  this  case  the  third  (SB)  arrival  reaches 
the  hydrophone  2/0  seconds  after  the  second  (B)  arrival  (see  Figure  3(b)).  From  this,  the  depth  of  the 
first  layer  is  given  by 


d\  =  jA/^vi 


where  A/2,3  is  the  time  difference  between  the  second  (B)  and  third  (SB)  arrivals,  and  v,  is  the  velocity  in 
layer  1. 

The  depth  of  the  second  layer  is  found  from  the  time  differences  between  the  direct  (D)  and 
bottom  reflected  (B)  paths.  If  the  angles  are  small,  then  tt  and  (t2- 13)  are  approximately  the  same,  and 
it  takes  the  bottom  reflected  path  2/3  seconds  longer  to  reach  the  hydrophone.  The  depth  of  the  second 
layer  is  then  given  by 


di  -  j  A'uv2 


where  A/1-2  is  the  time  different  between  the  direct  (D)  and  second  (B)  arrivals,  and  v2  is  the  velocity  in 
layer  2. 

In  general,  the  horizontal  distance  x  will  not  be  zero  (it  is  not  in  the  best  interest  of  the  ship’s 
crew  to  remain  over  the  site  where  an  explosive  charge  has  been  dropped).  This  method  will  give 
incorrect  depths  for  nonzero  distances  x  since  the  raypaths  are  not  vertical. 


Non-zero  ranges 

When  the  distance  x  is  not  zero,  it  is  necessary  to  correct  for  the  difference  between  the  vertical 
distances  and  d2  and  the  actual  length  of  the  raypaths.  Since  the  depths  dt  and  d2  are  unknown,  the 
angles  ot]  and  a2  are  also  unknown.  The  strategy  will  be  to  initially  use  the  values  d\  and  d2  found  by 
assuming  vertical  ray  paths.  Then,  given  x ,  corresponding  angles  cq  and  a2  may  be  found.  These  angles 
are  then  used  to  include  the  effect  of  non- vertical  ray  paths  in  the  calculation  of  depths  d ,  and  d2.  Since 
d ]  and  d2  are  not  exact,  oq  and  a2  are  not  exact,  and  the  calculated  dt  and  d2  will  not  be  correct,  but  will 
be  closer  to  the  correct  values  than  d\  and  d2.  This  procedure  may  be  iterated  using  the  calculated  values 
dj  and  d2  as  the  new  estimates  of  d\  and  d2  to  obtain  successively  closer  values  to  the  correct  depths. 
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Relations  for  finding  oq  and  a2  from  d\  and  d2  are  derived  from  the  geometry  of  the  ray  paths. 
From  Figure  4(a), 


tan(cq) 


x 

d\  +  ld2 


=  mUm<~d^2d ^ 


(1) 


From  Figure  4(b), 

*1 

tan(ot2)  =  — 
di 


and 


tan(ct2) 


l/2sa 
Wt  +  dz) 


(2) 


Equating  these  last  two  expressions  and  solving  for  x2  gives 

X\  \Hx2 

d\  (d  i  +  d£ 


2  (d\  +  d^  Xi 

d[ 


The  total  distance  is  x  =  x\+x2.  Using  this  to  findxj  gives 
Xi  =  x  -  x2 


Xi  =  x 


2  (d\  +  d  2)  X\ 


X  =  X\ 


1  + 


2(d\  +  d£ 

d[ 


X _ 

2(d,+d2) 


From  this  and  equation  (2), 


J,  ,  2W1  +  d2)l 


x 

(3  d\  +  2^2) 


tantaj 


Therefore, 


Estimated  depths  (assuming  vertical  raypaths)  d\  and  d2  may  be  substituted  into  (4)  to  obtain  estimates  of 
the  angles  oti  and  012.  These  angles  then  be  used  along  with  the  time  delays  to  find  depths  d\  and  d2- 
Since  these  depths  are  now  calculated  using  non-vertical  ray  paths,  they  will  be  closer  to  the  correct  shot 
and  seafloor  depths.  Expressions  for  d [  and  d2  in  terms  of  ,  a2 ,  and  the  time  delays  are  derived  next. 


From  Figure  5,  the  time  difference  between  the  third  and  second  arrivals  is 
A*2,3  =  (2t0+  2f3  + 14)  —  (2  / 1  -t-  f  2) 

/  2dx  2  d2  dx  \  /  2  d2  d ] 

vv!Cos(a2)  v2cos(02)  vjCosCotz)^  v'v2cos(a1)  v,cos(a!) 

( _ 3 _ 1  \  <*1  +  /  2  2  \ 

^  cosia^)  cos(a,K  v,  '  cos(ct2)  cos(a[)  '  v2 

aj  and  a2  are  known  (from  equations  (4a)  and  (4b)),  so  define  the  constants 

k  = _ 3 _ 1 

1  005(02)  cos(oq) 

=  -I _ 

cos(a2)  cos(«i) 


then  At  23  may  be  written  as 


The  difference  between  the  second  and  first  arrivals  is  (Figure  5) 


At  13  =  (2  f  i  +  1 2)  —  fa 


2  d2 

v2cos(cti) 


+ 


— rf-1— ) 
Vi  cos(oti) 


—  'Id2  +  x2 

Vl 


(1^1  +  - 1 }  1 
'  V2  Vj  cos(oti) 


±'Jd?T72 


Solving  (5)  for  d  t  gives 


di 


vi  ,  *2  Vi 

*7  M  "  *777 


(6) 


(7) 


The  expression  in  (7)  may  be  used  in  (6)  for  the  dx  enclosed  in  parentheses  (d,  under  the  radical  will  be 
replaced  later) 


At 


1 


U  - 


cos(cti) 


2^2  1  /  V!  k 2  v,  , 

-  +  — \~r~^23  ~  - d-2) 

v2  vj  kx  k 1  v2 


Vl 


V72 


1  +  *' 


1 


cos(ai) 


2  d2 
v2 


At 


2.3 


*2  d2 
*  1  v2 


-  —  ^ Id 2  +x2 


-  do 


v2cos(al)  v2cos(a1) ' 


-) 


At 


2.3 


^cosCa,) 


■<dT 


+  JT 


/  2  V] 

'  v2cos(a!) 


k2v. 


*1  v2cos(ct)) 


(  v  1  At 2,3 

'  /t]  cos(aO 


-  V!  At 


To  simplify  the  above  expression,  define 
2  vt  *2vi 

Ci  =  -  -  - 

v2cos(cti)  iiV2cos(a1) 


(8) 


k  i  cos(oq) 


and  (8)  may  be  written  as 


idl  +  x2  =  cxd2  +  c2 


Next,  substitute  the  expression  for  dx  from  (7)  into  (9)  and  square  both  sides. 

/  v,  k 2  vi  \2  ,  ,  ,  , 

A/2j  -  -r-—d2)  +  x  =  ci  d2  +  2c!C2d2  +  c2 

*1  <i  v2 


.viA/2.3.2  2v2  k2At2t3  /t2Vn2.2  2  2  ,2  „  2 

(  r  )  -  71 - ^2  +  (t — ~)  d2  +  x  -  ci  d2  +  2 c\c2d2  +  c2 

k i  t  v2  *iv2 


tk2v\\  ^  ,?  2v2  k2Atv  ( V] Af23  ^2 

(7— )  -  Cl  ^2  + - - 2c,c2  d2  +  - )  -ci+X2\  = 


This  is  in  the  form  a  d2  +  b  d2  +  c  =  0  and  the  quadratic  formula  may  be  used  to  solve  for  d2. 
Once  d2  is  known,  d !  may  be  found  from  equation  (7) 

vi  k2  v, 

d'  =  ~  T~~Td2 

*1  *1  v2 


This  estimate  of  the  shot  depth  (d,)  and  seafloor  depth  (di  +  d?)  includes  the  effect  of  non¬ 
vertical  ray  paths  since  angles  oq  and  a 2  were  used  in  the  calculation.  These  values  of  dx  and  d2  are  then 
used  to  calculate  new  values  for  a!  and  o^,  and  the  process  is  repeated.  dx  and  d2  should  converge  to 
constants  within  a  few  iterations. 
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3.  Algorithm 


Required  inputs 

a)  differences  in  time  of  arrivals  At x2  and  ^2j- 

b)  distance  the  ship  has  traveled  from  the  drop  site  at  the  time  of  detonation. 

x  =  (ship’s  speed)  x  (detonation  time  -  time  charge  was  launched ). 

c)  average  sound  speed  V!  and  v2  in  layers  1  and  2. 


Outputs 

a)  shot  depth  (d  0 

b)  distance  from  shot  to  sea  floor  ( d 2). 
Bottom  depth  is  then  given  by  dt  +  d2- 


1.  Find  approximate  values  for  dx  and  d2  by  assuming  vertical  ray  paths 

j  1  . 

d\  =  TA,13V1 


d2  ~  YA,1.2v2 


2.  Calculate  a(  and  a2 


ai  =  arctan(— — - — —  ) 
dx  +  2d2 


a2  =  arctan( — — — ) 
3  dx  +  2d2 


jS 


& 


3.  Calculate  constants  k  i  and  k2 


*1 


3 _ 1 

cos(a2)  cos(a,) 


_ 2 _ 2_ 

COS(Ct2)  COS(a!) 


4.  Calculate  constants  c !  and  c2 

2vj  k2vx 

C  =  -  —  - 

v2cos(a,)  ^vjCosCcq) 


C2 


vl^f2.3 

jkjCOsCaj) 


vi 


5.  Calculate  coefficients  for  the  quadratic  formula 


/*2yl\  2 


b  =  - 


2  vf  k2  Ar2-3 
*i2v2 


-  2ci c2 


/  v  1  A/^3  ^  2 

c  =  )  -  c|  +  JSZ 


6.  Find  depths  d  2  and  d  x 

-b  ±  -  4ac 


d2  = 


2a 


,  V1  a  *2  Vi 

dl  -  I7a,w  "  *7^2 


7.  Use  d !  and  d2  found  in  step  6  as  the  next  dx  and  d"2  in  step  2  and  repeat. 

After  dx  and  d2  have  converged  to  constants,  dx  gives  the  shot  depth  and  (d,  +  dO  gives  the  sea  floor 


4.  Example 

A  Fortran-77  program  which  implements  the  algorithm  on  page  12  is  given  in  the  Appendix. 
To  illustrate  its  use,  shot  and  seafloor  depths  will  be  found  using  the  data  presented  in  Figure  2.  This  par¬ 
ticular  data  was  obtained  from  a  refraction  experiment  conducted  on  the  Monterey  Deep-Sea  Fan  in  1981 
by  the  Marine  Physical  Laboratory.  The  research  vessel  FLIP  (FLoating  Instrument  Platform)  was 
moored  in  2800  meters  of  water  and  utilized  as  a  receiving  station.  A  20  element  hydrophone  array  was 
deployed  at  mid-depth  in  the  water  column.  The  output  of  nine  of  these  hydrophones  along  with  the 
hydrophone  streamer  data  is  shown  in  Figure  2. 

A  Navy  tug  steamed  away  from  FLIP  and  dropped  explosive  charges  (Mk  94  SUS  +  256# 
TNT)  into  the  water  at  ranges  between  3.5  km  to  37  km  from  FLIP.  As  each  shot  was  launched  the  radar 
range  between  FLIP  and  the  tug,  and  the  time,  were  recorded.  Velocity  of  the  tug  was  obtained  from  this 
information.  The  distance  traveled  from  the  drop  site  by  the  tug  when  the  detonation  arrived  at  the  strea¬ 
mer  hydrophone  could  then  be  calculated. 

Mean  velocities  were  obtained  by  integrating  the  velocity  profile  between  the  surface  and  1820 
meters  for  layer  I,  and  from  1820  meters  to  2800  meters  for  layer  2. 
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The  following  input  to  the  program 


Enter  mean  velocity  in  layer  1  (meters/sec) 

1485 

Enter  mean  velocity  in  layer  2  (meters/sec) 

1492 

Enter  first  and  second  arrival  time  difference  (sec) 
1.19 

Enter  second  and  third  arrival  time  difference  (sec) 
2.18 

Enter  distance  ship  traveled  from  drop  site  (meters) 
2738 

Enter  maximum  number  of  iterations  to  perform 
10 


produces  the  output 


shot  depth:  1788.4 
sea  floor  depth:  3016.6 


These  are  reasonable  values  for  this  experiment.  The  shot  depth  of  1788  meters  is  within  the 
rather  large  uncertainty  of  the  detonation  device.  A  bathymetric  map  of  the  area  shows  the  sea  floor 
depth  to  be  between  2800  meters  at  FLIP  and  3200  meters  at  a  range  of  32  km  from  FLIP. 
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Appendix:  Fortran-77  program 
Implements  the  algorithm  on  page  12 
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c  - - 

c  Iterative  method  to  find  shot  depth  (dl)  and  seafloor  depth  (dl  +  d2) 
c  from  hydrophone  streamer  data. 

c  - 

real  kl,  k2 


c  - Input - 

write(6,*)  ’Enter  mean  velocity  in  layer  1  (meters/sec)’ 
reaa(5,*)  vl 

write(6,*)  ’Enter  mean  velocity  in  layer  2  (meters/sec)’ 
read(5,*)  v2 

write(6,*)  'Enter  first  and  second  arrival  time  difference  (sec)’ 
read(5,*)  dtl2 

write(6,*)  ’Enter  second  and  third  arrival  time  difference  (sec)’ 
read(5,*)  dt23 

write(6,*)  'Enter  distance  ship  traveled  from  drop  site  (meters)’ 
read(5,*)  x 

write(6,*)  'Enter  maximum  number  of  iterations  to  perform’ 
read(5,*)  iter 


c  Step  1. 

c  Estimate  dl  and  d2  by  assuming  ray  paths  are  vertical 
dhatl  =  .5  *  dt23  *  vl 
dhat2  =  .5  *  dt!2  *  v2 


a  =  (k2*vl/(kl*v2))*(k2*vl/(kl*v2))  -  cl*cl 
b  =  -2.* v  1  * v  1  *k2*dt23/(k  1  *k  1  * v2)  -  2.*cl*c2 
c  =  x*x  -  c2*c2  +  (vl*dt23/kl)*(vl*dt23/kl) 


Step  6. 

Note:  The  solution  to  the  quadratic  formula  is  -b  +/-  sqrt ... 
For  the  data  used  in  this  example  the  minus  sign  gives 
the  correct  result 

d2  =  ( -b  -  sqrt(  b*b-4.*a*c  ) )  /  (2.*a) 
dl  =  vl*dt23/kl  -  k2*vl*d2/(kl*v2) 


Make  these  values  of  dl  and  d2  the  new  estimates,  and  repeat 
dhatl  =  dl 
dhat2  =  d2 

100  continue 


- —  Output . 

write(6,*)  '^not  depth:  ’,dl 
write(6,*)  'sea  floor  depth:  dl+d2 


